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[l] A perennial storage of water in a fim aquifer was 
discovered in southeast Greenland in 2011. We present the 
first in situ measurements of the aquifer, including densities 
and temperatures. Water was present at depths between ~12 
and 37 m and amounted to 18.7 ±0.9 kg in the extracted 
core. The water filled the fim to capacity at -35 m. 
Measurements show the aquifer temperature remained at the 
melting point, representing a large heat reservoir within the 
fim. Using model results of liquid water extent and aquifer 
surface depth from radar measurements, we extend our in 
situ measurements to the Greenland ice sheet. The estimated 
water volume is 140 ±20 Gt, representing ~0.4 mm of sea 
level rise (SLR). It is unknown if the aquifer temporary 
buffers SLR or contributes to SLR through drainage and/or 
ice dynamics. Citation: Koenig, L. S., C. Miege, R. R. Forster, 
and L. Bmcker (2014), Initial in situ measurements of perennial 
meltwater storage in the Greenland fim aquifer, Geophys. Res. Lett., 
41, doi: 1 0. 1002/20 1 3GL058083. 

1. Introduction 

[ 2 ] Mass loss from the Greenland ice sheet (GrIS) is acceler- 
ating at a rapid pace [Shepherd etal., 2012; Hanna etal., 2013]. 
In the last decade, surface mass balance processes, dominated 
by meltwater from the ice sheet, account for over half of the 
mass loss [van den Broeke et al., 2009]. Very little, however, 
is known about the quantity of meltwater that flows through 
the hydrologic system, its pathways, or the rates of refreezing 
and retention that occur as the meltwater travels toward the 
sea [Harper et al., 2012; Humphrey et al., 2012; van Angelen 
et al., 2013], Studies of meltwater retention and refreezing in 
fim are limited in number and location mostly to Western 
Greenland [Boggild et al., 2005; Humphrey’ et al., 2012; 
Harper et al., 2012; Rennermalm et al., 2013] though simula- 
tions show that 42 ± 4% of rain and meltwater is either retained 
or refrozen in the fim without reaching the sea [ van Angelen 
et al., 2013]. Only recently, meltwater was observed to be stored 
in the fim over the winter in a large (70 ±10 x 10 3 km 2 ) fim 
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aquifer, concentrated in Southeast Greenland where both accu- 
mulation and melt rates are high [Burgess et al., 2010; Forster 
et al., 2013]. The fim aquifer is a perennial storage of water at 
depth, between 5 and 50 m below the surface, controlled by local 
topography, seasonal snowpack thickness, and melting and 
refreezing rates. Fim aquifer analogs are found in temperate gla- 
ciers [e.g., Fountain, 1989] and ice caps [e.g., de Woul et al., 
2006], but the Greenland aquifer is deeper, larger, and subsists 
year after year [Forster et al., 2013]. Here we present the first 
density and temperature profiles associated with the 
Greenland fim aquifer; we estimate the volume of water stored 
and discuss its implications for sea level rise (SLR) in terms of 
the ice sheet mass budget and water retention. 

2. Field Measurements and Analysis 

[ 3 ] During the first 10 days of April 2013, we took in situ 
measurements of the aquifer at a site, FA 13 (66.1812°N, 
39.0435°W, 1563 m above sea level (asl)), located between 
the two Arctic Circle Traverse (ACT) core sites (ACT1 1-A 
and ACT1 1-A2) where the water-saturated fim was first ob- 
served (Figure 1) [Forster et al., 2013]. The ACT traverses, 
conducted in the spring of 20 1 0 and 2011, measured accumu- 
lation rates in the high-accumulation region of SE Greenland 
using a 400 MHz ground-penetrating radar and fim cores 
[Miege et al., 2013]. The closest seasonally dry fim core to 
our site is ACT11-B (66.2174°N, 39.5658°W, 1806masl), 
located -24 km West of our site and -7.5 km West of the lateral 
edge of the fim aquifer. ACT11-B, referred to as the dry fim 
core, has well-defined melt layers throughout the 60 m core 
comprising -17% of the core length. At site FA13, the water 
layer was mapped by radar in 2011 to be close to the surface 
at a depth of 15 m. Two cores, FA13-A and FA13-B, located 
2 m apart, were drilled from the surface to depths of 60.4 m 
and 32.0 m, respectively, into the aquifer to measure water 
level, density, and temperature. A mechanical drill was used 
in the seasonally dry fim above the water table (upper -10 m), 
and an electrothermal drill was used to drill through the 
saturated fim. 

2.1. Density Measurements 

[ 4 ] Figure 2a shows the gravimetric density profiles from 
cores FA13-A and FA13-B at -20 cm vertical resolution 
and from the closest seasonally dry firn core, at -1 m vertical 
resolution. (Densities were measured to a depth of 50.5 m in 
FA13-A.) Measurement errors of±5% are assumed due to 
core volume irregularities. Additional unquantifiable density 
errors are introduced when water drains from the saturated 
firn cores during extraction and handling, giving a low bias 
to all densities with water saturation (see supporting informa- 
tion). Density profiles of the seasonally dry fim above the 
water table (top -12 m) mimic that of the closest dry fim core 
and clearly show higher-density layers characteristic of 
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Figure 1. Map of 2013 study area. Cores FA13-A and FA13-B (blue circles) are located between ACT1 1-A and ACT1 1-A2 
(green circles) where the water table was first observed. ACT1 1-B (green square) is the closest dry firn core. Elevation from 
Howat et al. [2013], 


infiltrated water from summer melt or rain and subsequent 
refreezing. Starting at ~8 m below the surface, the seasonally 
dry fim has a very well-connected percolation drainage 
network evident by piping structures seen within the core 
(Figure SI in the supporting information). Water was first 
encountered at a depth of 12.2±0.1m in both boreholes. 
Water immediately filled the boreholes (~12 cm diameter) be- 
low 12.2 m (Figure S2) and remained at this level over the 
4 day sampling period. When the second borehole (FA13-B) 
was drilled, the water level did not change in FA13-A, imply- 
ing that the firn aquifer is unconfined [Fountain, 1989] and 
gravitationally fed by surface meltwater. 

[ 5 ] Core segments between 12.2 m and 36.9 m were com- 
posed of a mixture of well-bonded fim and well-connected 
pore space either filled or partially filled with water evident 
by the densities above that of ice (Figure 2). FA13-A and 
FA13-B density profiles were nearly identical and were com- 
bined into a single-FA13 density profile for further analysis 
(see Figure 2 and supporting information for details). Due 
to the lack of an adequate depth densification model for the 
percolation zone with water-saturated fim, two fits were 
used in our analysis. The measured densities were fitted with 
an exponential curve introduced by Horhold et al. [2011] 
(Figure 2a) and a modified Herron and Langway (HL) fit 
[Herron and Langway, 1980] (Figure 2b). The HL fit is 
empirically based for dry snow compaction; however, here 
only the form of the equation is maintained. The critical ice 
density, point at which recrystallization and deformation 
are the dominant densification processes, was allowed to 
vary. Therefore, in this melt-/freeze-dominated environment 
the model parameters based on dry snow processes have no 
physical meaning (i.e., the best fit mean annual temperature 
and accumulation are not reasonable for this location on 
the ice sheet). Both fits were adjusted by minimizing the 
root-mean- square error (RMSE) and were nearly identical 
in fit (exponential fit, RMSE = 43.7 kg/m 3 and HL fit, 
RMSE = 44.3 kg/m 3 ). The exponential fit was used in our 


analysis because it includes an additional exponential term 
over HL fit, thus allowing more flexibility in representing 
the measurements. 

[6] While we had no direct way of quantitatively determining 
pore close-off depth in the field, we used the simple techniques 
of manually blowing air into the 20 cm core sections, visual 
inspection, and visible photography. We estimated a pore 
close-off depth between 33 and 35 m (Figure 2a). In dry fim 
conditions pore close-off is defined as ~830kg/m 3 [Cuffev 
and Paterson, 2010], which is the density at which air bubbles 
become entrapped in the ice. However, with water-saturated 
fim, where densities are often over 900 kg/m ' (Figures 2 and 
S3) [Fountain, 1989], a new metric is needed to define the 
lower boundary of the fim aquifer. We define the pore water 
refreeze depth as the depth where water in the pore space 
produces a maximum density; below this depth the water in 
the pore space refreezes, causing a gradual decrease in density 
to the density of ice. The pore water refreeze depth occurs at 
36.9 m detennined by the exponential fit (Figure 2a). At this 
depth the maxima-fitted density (9 15 kg/m 3 ) is reached 
and densities below this depth show that the small amounts 
of entrapped water are refreezing to ice. We define the 
bottom of the fim aquifer as the pore water refreeze depth 
(rounded to 37 m). Though it is certain that small amounts 
of water exist below our defined bottom of the fim aquifer, 
it is refreezing to ice. 

2.2. Temperature Measurements 

[ 7 ] Thermistor strings with 0.5 m and 2.5 m spacing between 
sensors were inserted in the FA13-B and FA13-A boreholes, 
respectively; each borehole was backfilled using fine surface 
snow, and the thermistors were left to record subdaily tem- 
perature profiles for a year or longer. Figure 3 shows the 
0.5 m temperature profiles from the FA13-B borehole. The 
surface of the water table is easily distinguishable at ~12m 
with a constant 0.1±0.1°C temperature recorded from this 
point to the bottom of the thermistor string. The temperatures 
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Figure 2. (a) Density profiles from ACT1 1-B, FA13-A, and FA13-B with exponential fits [Horhold etal., 2011]; gray area is 
the field-observed pore close-off depth, and horizontal black line is the surface of the water table, (b) Density profiles from 
ACT11-B and combined FA13 with a modified Herron and Langway [1980] fit. The horizontal lines show the surface of the 
water table and the pore water refreeze depth. The black vertical lines from left to right represent the infiltration ice density 
[Harper et ah, 2012] and ice density. 


measured are indistinguishable from the melting point and 
assumed 0°C with both water and ice present. In contrast, tem- 
peratures recorded in the seasonally dry fim above the water 
table show surface-forced warming in the spring/summer. 
The first surface-melting events were recorded on 10 June 
2013. The temperature profiles highlight the temperate 
water/ice conditions at depth, much warmer than expected 
for the region. The aquifer, through water percolation and 
phase changes, represents a heat reservoir in the GrIS which 
is not currently accounted for in fim models. 

2.3. Annual Stability of the Water Table Elevation 

[s] Annual variations in the depth of the water table sur- 
face were measured from 2011 to 2013 using ground-based 
and airborne radars, center frequency of 400MFIz and 
750MFIz, respectively (Table 1). For this assessment we 
use only the depths calculated from radar measurements 
and ice sheet surface elevations. We neglect change from 
variations in fim compaction and ice flow. The accuracy of 
the radar-derived depths is ± 0.3 m based on the vertical reso- 
lution of the radars. These data show a constant surface 


elevation of the water table at this location from 2011 to 
2012. In 2013, the water surface elevation rose by 0.7±0.3 m, 
coinciding with a surface lowering of 2 m. This occurs after 
the 2012 extreme summer melt over the GrIS [Nghiem et al., 
2012; Hall etal., 2013]. While the time series is short, 3 years, 
the annual variation of the water table surface, as summarized 
in Table 1, shows that between April 2012 and April 2013 the 
aquifer likely filled caused by the intense melt season and 
identified by the rise in water table. 


Table 1 . Annual Variations of Water Table Surface 


Radar Data 

Collection 

Date 

Surface 
Elevation (m) 

Water Table Depth 
From Surface (m) 

Water Table 
Surface Elevation (m) 

2 May 2011 

1565.5 

15.0 

1550.5 

17 April 2012 

1565.6 s 

15. l b 

1550.5 

3 April 2013 

1563.5 

12.3 

1551.2 


a From Airborne Topographic Mapper data [ Krabill , 2012], 
“From Accumulation Radar data [Leuschen, 2012], 
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Figure 3. Daily averaged temperature profiles from FA13-B 
borehole. The horizontal black line represents the surface of 
the water table. Each thermistor is represented by a black dot. 

3. Water Saturation and Implications for Sea 
Level Rise 

[ 9 ] To approximate the volume of stored water from the 
new cores, it is assumed that the pore space available for 
water saturation is equal to the pore space calculated from 
the closest dry firn density profile and that the density 
difference between the water-saturated core and the dry 
core is attributed to water (see supporting information). 
Figure 4 shows the percentage of pore space available 
for water saturation and the percentage of water at this 
site. At 33.5 m the percentage of water and percentage of 
available pore space overlap, showing complete saturation 
of the fim, ~3.5m above the pore refreeze depth. The 
field-estimated pore close-off depth (33-35 m) overlaps 
with water saturation where compaction and water had likely 
decreased the ability for air to move through the core. 

[ 10 ] The water table is 24.7 m thick at this location, and the 
core contained 1 8.7 ± 0.9 kg of water, from its water-saturated 
surface to the pore water refreeze depth. The fim aquifer 
surface mapping done by Forster et al. [2013] showed a 
Gaussian distribution of depths with a mean of 23 m. 
Assuming that the bottom of the entire aquifer is constant 
and stable at 37 m depth, this gives a mean thickness of 
14 m with a two sigma thickness range of 0 to 28 m. We cal- 
culate the total volume of the water-saturated aquifer, using 
the mean thickness (14 m) and the modeled area over the 
GrIS ice sheet (70±10 x 10 3 km 2 ) [Forster et al., 2013], 
to be 980 ± 140 km 3 . Assuming the same density of water 
measured at the core site throughout the entire firn aquifer, 
the total estimated amount of water stored in the aquifer 


throughout the ice sheet is 140 ±20 Gt, which represents 
~0.4 mm of SLR. 

[ 11 ] It is unknown from our measurements if the stored 
water is connected to the ice sheet basal hydrology or if it 
has any pathways to drain into the ocean. Two hypotheses 
present possible pathways for this stored water to exit the 
aquifer. (1) The aquifer is connected to a well-developed 
englacial hydrologic network, including crevasses and 
moulins that drain some portion of the aquifer at a relatively 
constant rate (seasonally) to the bed, similar to surface melt 
in Western Greenland [e.g., Zwally et al., 2002]. This hy- 
pothesis assumes that the water table has persisted over a 
long enough time scale to develop a hydrologic network. 
Additional studies on water residence time, flow paths and 
rates, and volume input are required to favor this hypothesis 
as well as a connection to crevasses in the region which has 
not yet been found. Or (2) the aquifer is primarily storing 
water in available firn pore space and will not release water 
until the pore space is completely saturated and/or a thres- 
hold is met, leading to a catastrophic release event, similar 
to jokulhlaups [e.g., Nye, 1976]. The second hypothesis is 
supported by our initial measurements showing available 
pore space (Figure 4) and a rise in the water table after 
record surface melting (Table 1). These hypotheses are not 
independent, and both scenarios need additional study. 
The firn aquifer likely plays a significant role in SLR by 
either accelerating glacial flow in the region and/or by 
storing and then catastrophically supplying large amounts 
of fresh water if connected to the ocean. 

4. Firn Capacity Within Firn Aquifer Regions 

[ 12 ] Assuming the second hypothesis, the firn aquifer 
may serve as a long-term reservoir for meltwater that tem- 
porarily buffers SLR. Flarper et al.[ 2012] used firn capacity 
(in kg/m 2 ) to quantify the capacity of the GrIS to buffer SLR 
by storing percolated meltwater in the firn that refroze. 
While the calculations done by Harper et al. [2012] are 
valid for most of the percolation zone, it is clear that in 
regions of high accumulation and high melt rates, like those 
in southeast Greenland where the PFA has been mapped, 



Water Filled Pore Space (FA13) 

Pore Space (ACT 1 1 B) 


_ 60 l i i i * * 1 1 j 

-10 0 10 20 30 40 50 60 70 

Percent 

Figure 4. Percent of pore space (ACT 11 -B) and percent of 
water (FA13) assuming equivalent fim porosity between the 
dry and water-filled cores. 
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Harper et al.’s [2012] capacity calculation must be adjusted 
for water retention which adds capacity. After adjusting 
calculations in aquifer regions, water storage represents an 
increase of 8.9% over regions of refrozen meltwater due to 
the difference in density between water and ice (see supporting 
information for details and Figure S4). Fim aquifer regions, 
therefore, allow for additional storage of meltwater, and there- 
fore a larger SLR buffer than accounted for by Harper et al. 
[2012], Though the fim aquifer can store more meltwater, 
the site studied is already filled to 50% of its capacity com- 
pared to the nearby seasonally dry location (see supporting 
information and Figure S5). Fim capacity is reached at 
~35 m showing that the bottom of the aquifer has reached full 
capacity. Future meltwater must continue filling the upper 
portion of the aquifer until either the fim is full, no longer 
allowing percolation of surface melt, forcing melt to flow 
directly to the ocean or a tipping point is reached where the 
stored water catastrophically releases [e.g., Nye, 1976], 

5. Conclusion 

[13] The perennial fim aquifer is large gravitationally fed 
liquid water body in the GrIS. At our study site the top of 
the water table is at 12.2 m and water extends for 24.7 m to 
the pore water refreeze depth at 36.9 m, though small 
amounts of water do exist below this depth. A longer time se- 
ries is needed to confidently correlate depth variations in the 
water table surface with specific melt and accumulation 
events, though there is evidence of a rise in the water table af- 
ter the extreme summer melt of 20 1 2. It is unknown if the fim 
aquifer is currently buffering SLR, with an additional 8.9% 
storage over other percolation regions, or if the fim aquifer 
episodically accelerates SLR by influencing ice dynamics 
or periodically draining to the ocean. 

[14] Indeed, quantitative measurements of water saturation, 
flow velocity, and permeability are needed to better character- 
ize the water volume at several additional sites across the fim 
aquifer. Further studies and modeling of the fim aquifer are 
also needed to fully understand its impact on SLR and to 
improve future predictions. These initial measurements do 
show that the volume of liquid water storage in SE 
Greenland is large. With an estimated ~0.4 mm contribution 
to SLR contained in the fim aquifer, it is an important source 
of previously neglected heat and water storage in the GrIS. 
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